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Abstract 

We examine the relation between the damping rate of a chiral fermion mode 
propagating in a hot plasma and the rate at which the mode approaches equilibrium. 
We show how these two quantities, obtained from the imaginary part of the fermion 
self-energy, are equal when the reaction rate is defined using the appropriate wave 
function of the mode in the medium. As an application, we compute the production 
rate of hard axions by Compton-like scattering processes in a hot QED plasma 
starting from both, the axion self-energy and the electron self-energy. We show 
that the latter rate coincides with the former only when this is computed using the 
corresponding medium spinor modes. 

PACS numbers: ll.10.Wx, 11.15.-q 
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1 Introduction 



As is well known, the properties of a chiral fermion propagating in a medium are modified 
from those in vacuum. Interactions give rise to the appearance of collective modes char- 
acterized by dispersion relations and damping rates. To be precise, if = (u,p) is the 
four-momentum of a fermion mode in the rest frame of the medium, the corresponding 
pole of the propagator is located, in general, at a complex value 

w = u p -i^, (1) 

where u p and 7 are real functions of p and correspond to the energy and damping rate 
of the mode, respectively. In particular, m/ = ui p =o can be interpreted as the mass of 
the fermionic excitation and 7 _1 , in accordance with linear response theory gives the 
characteristic time scale for the decay of the mode. In a field theoretical description, u p 
and 7 are determined (for weak damping) from the real and imaginary part of the fermion 
self-energy, respectively. 

We recall, however, that there exists a distinct concept also related to the imaginary 
part of the fermion self-energy. This quantity, which we will call the total reaction rate T, 
is interpreted in such a way that its inverse gives the time scale for the fermion distribution 
to approach equilibrium 0, |!|. It would then appear that one has two different physical 
quantities, the damping rate 7 and the total reaction rate T, both obtained from the 
imaginary part of the fermion self-energy. This has indeed been a source of some confusion 
in the literature, for instance, in a paper on the subject, two different expressions are used 
to compute the damping and total reaction rate in terms of the self-energy [|J . 

The equivalence between the above two quantities was established by D'Olivo and 
Nieves for a chiral fermion interacting with a scalar and a massive fermion through a 
Yukawa coupling. They showed that in this case 

r = 7, (2) 

provided T is defined using the appropriate wave functions of the chiral fermion mode 
in the medium. In this paper we show explicitly how this equivalence is realized in the 
realm of hot gauge theories where the resummation method of Braaten and Pisarski || has 
provided a consistent framework for the computation of the leading-order contributions 
to damping rates. We also stress the role played by wave function renormalization of the 
medium spinor modes for the above equivalence. We illustrate this role by computing the 
production rate of hard axions in a hot QED plasma through Compton-like scattering 
processes, showing that the rate obtained by considering the axion self-energy coincides 
with the rate obtained from the expression for the electron self-energy only when the 
appropriate medium spinors are used. 

This paper is organized as follows: Sections II and III are devoted to a brief recollection 
of some well established concepts regarding the self-energy of a chiral fermion in a medium 
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and its relation to the fermion damping rate. In section II we review how the damping rate 
can be written in terms of the imaginary part of the self-energy and the spinors satisfying 
the effective Dirac equation in the medium. In section III we briefly recall how a complete 
leading order calculation yields gauge independent poles for the fermion propagator in the 
medium in a covariant gauge. In section IV we express the imaginary part of the self- 
energy as a total reaction rate. We identify the physical processes involved and establish 
the equivalence between the damping rate and the total reaction rate when the latter is 
expressed in terms of the spinors in the medium. In section V we compute the production 
rate of hard axions in a hot QED plasma. We finally discuss our results in section VI. A 
short appendix collects notations and definitions of some quantities appearing in section 
IV. 

2 The damping rate 

In this section we express the damping rate in terms of the self-energy, following D'Olivo 
and Nieves ||, thereby establishing our notation and preparing for the discussion in the 
next section. The inverse propagator of a chiral fermion in a medium can be written in 
general as 

^- 1 (P) = ^-S, (3) 

where £ is the effective self-energy induced by the medium. Hereafter, capital letters are 
used to refer to momentum four- vectors. At T = the theory is chirally invariant, due to 
the vanishing of the fermion mass. Chiral invariance is unaffected by the presence of the 
medium, hence the inverse propagator in the medium rest frame can be expressed as 

iS-\P) = A olo -A s ^-p 

= Ia; 1 (p)( 7o + 7-p) + ^a: 1 (p)(7o-7-p) (4) 

in terms of the functions A±(P) = (Aq =p A s ) _1 , where we have denoted the unit vector 
p/p by p- The propagator poles are given by A^ X (P) = 0. In the positive-energy sector, 
Eq. (||) yields four propagating modes, two with positive and two with negative helicity 
over chirality ratio 0. For definitiveness, let us work with a negative-helicity solution 
and look for the positive-energy poles in A + (P). The same arguments apply to negative- 
energy and positive- helicity solutions and the poles of A_(P). 

A^ 1 (cj,p) has in general a real and an imaginary part. Writing uj as in Eq. (fj), the 
equation determining the poles is 

ReA^ 1 ^ - ij/2,p) +«ImA; 1 (w p - ij/2,p) =0. (5) 

In the following, we consider the case of weak damping, 7 <C u p , where the physical 
concept of a propagating mode is still meaningful. Eq. (|5|) can then be solved by expanding 



3 



it in powers of 7 and retaining terms at most linear in 7 and Im A + x . At zeroth order, u p 
is obtained from 

ReA-\u p ,p) = 0, (6) 
whereas 7 is given by the first-order equation, which yields 



7 

where 



- = Z p ImA + 1 (^p,p) , (7) 



-1 



The factor Z p coincides with the residue at the pole of the one-particle contribution to 
the propagator and it is thus the normalization factor that has to be taken into account 
to construct the spinors representing one-particle states ||. In the context of many-body 
physics, Z p represents the probability factor that needs to be considered for the various 
processes involving the fermion mode M. 

In order to relate 7 to ImS, recall that the one-particle states are represented by 
spinors satisfying the effective Dirac equation 

(P-ReE)u(P) = (9) 

obtained by neglecting the absorptive part of the effective self-energy. Real and imaginary 
part of £ are defined, as usual, by 

Re £ = i(£ + 7 £t 7 
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ImE = ^ (S -7°S t 7°) • (10) 

In the following, we will not need the explicit solutions of Eq. but rather work with 
the corresponding projection operators. Choosing a normalization corresponding to that 
of the one particle states M, the projector for the solution under consideration can be 
written as 

u(P)u(P) = Z p iu p Lfl, (11) 

where L — (1 — 7s)/2 and is a four- vector with components = (l,p) in the plasma 
rest frame. From Eqs. ([]), (|]) and (0) we get, by a straightforward computation, 

7 = — ^Tr [/HmS^p)] 

= -—u(P)lm^(co p ,p)u(P), (12) 
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where in the last line we have used Eq. 
evaluated at uj 
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p- 



Eq. (12 



O). Note that in Eq. flT2| ) the self-energy is 
is valid independently of any approximation scheme, 
perturbative or not, and it only requires that the damping rate be small compared to the 
energy of the mode. This equation was already used in Ref. || to establish the equality 
of damping and reaction rates in a model with Yukawa coupling. In the next section, we 
review how the damping rate is obtained in hot gauge theories. 



3 The HTL self-energy 

In a gauge theory, a plasma of massless particles is characterized by two quantities, the 
temperature T and the dimensionless coupling constant g. If g <C 1 there exist two dis- 
tinct energy scales, T and gT. The first is characteristic of individual particles, whereas 
the second is typical of collective phenomena. In the Hard Thermal Loop (HTL) approx- 
imation, bare perturbation theory is used to compute the leading contributions to the 
self-energies and vertices. For soft external momenta (i.e. of order ~ gT), these are of 
the same order of magnitude as the bare quantities that they modify, thus the necessity 
of resummation. HTLs provide effective propagators and vertices necessary for a consis- 
tent perturbative expansion. In the following, we will consider exclusively soft external 
momenta. The lowest order HTL contribution to the fermion propagator can be shown 
to be real (for time- like momenta) and gauge-independent ||. A complete leading order 
calculation of the imaginary part of the self-energy requires, however, the sum of the two 
diagrams depicted in Fig. 1. The heavy dots indicate effective propagators and vertices. 
Fig. la is the usual diagram for the fermion self-energy modified by the use of effective 
propagators and the effective two-fermion one-gauge boson vertex. Fig. lb is an addi- 
tional contribution coming from the effective two-fermion two-gauge boson vertex, upon 
contracting the gauge boson lines. 

Let us denote by *E(P) the effective self-energy calculated to leading order in its real 




(a) (b) 

Figure 1: The self-energy graphs for fermions at one-loop level in the effective expansion. 
Heavy dots denote the effective propagators and vertices arising in the HTL approximation 
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and imaginary parts. In the imaginary time formalism of Thermal Field Theory, the 
contributions from Figs, la and lb are, respectively 

*£ a (P) = -g^C F TY,J-0^r,(P,P-K)D, u (K)S(P-K)r u (P-K,P) 

*E*(P) = -jC F T^J-^r^(P,P,K)D, v (K), (13) 

where Cp is the Casimir invariant corresponding to the fundamental representation of 
the gauge group. The effective two-fermion one-gauge boson and two-fermion two-gauge 
boson vertices are denoted by Y ^ and T^ u , respectively. D^ u and S represent the gauge 
boson and fermion effective propagators, respectively. 

In a covariant gauge, D^ u has a piece proportional to the gauge parameter £. Let us 
separate its contribution to the effective self-energy writing 

*£(P) = *£ (P) + *£ 5 (P) (14) 

where *£ (P) is gauge-parameter independent and *E^(P) is the contribution to Eq. (|T3| ) 
from the term proportional to the gauge parameter in D^ u , namely, from 

e k,k v 

K 2 K 2 ' 

By using the Ward identities satisfied by the HTLs we can express the gauge-parameter 
dependent term in Eq. (|T4|) as 

^{P) = a{P)S-\P), (15) 

where we have defined 



a{P)^ig 2 C F TY. I 



d 3 k 



S(P - ^S-^P) 



(K 2 ) 2 (K 



2\2 



(16) 



The function a(P) is infrared singular. In what follows we assume that a proper infrared 
regulator has been introduced and that this will be taken to zero only until the very end 
of any calculation flTU |. 

The modification that *£(P) introduces in the fermion propagator can be written, 
after analytical continuation to Minkowski space, as in Eq. (H). Using Eqs. (0) and fllBD , 
we have 

i *S~ 1 (P) = P- *£ (P) - i a{P)S-\P) . (17) 

A simple power counting shows that a(P) is of order g. Therefore, to the given order, it 
is allowed to replace S' _1 (P) by *S' _1 (P) in Eq. ( |17| ) which now becomes 

t*S-\P) = P-*£(P) 

[l-a(P)][P-^ (P)], (18) 
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whence 



*E(P) = a{P)P + [1 - a(P)} *£ (P) . (19) 

Notice that in principle, there are other non-leading contributions to the real part of 
*£(P) which are not taken into account in Eq. fll9|) . In fact, for the calculation of the 
damping rate to leading order, it is only necessary to consider the leading contributions 
to the real part of *E(P). Nevertheless, in the following we keep a non-leading term 
containing a(P) to illustrate that its presence does not affect the gauge independence of 
the damping rate. 

Eq. (|T8|) implies that the poles of the propagator are gauge-parameter independent, 
since a pole of [P— *E (P)] _1 is also a pole of *S. Within our approximation, the effective 
Dirac equation for the one-particle fermion mode is 

[P-Re*£(P)]*<u( p ) = 0. (20) 

By neglecting the absorptive with respect to the dispersive part of *£ (P), Eq. (|20|) 
becomes 



[P-Re*E (P)]XP) = 0. (21) 

According to Eq. (|TT|) , a factor [1 — Rea(P)]~ 1//2 has to be absorbed into the normalization 
of the spinor *u(P) since this factor is part of the residue of the fermion propagator at 
the pole. 

From. Eqs. (O) and ([HJ), the damping rate 7 can be written as 
7 = — *u(P)Im*£(cj p ,p)*u(P) 

U p 

= —[l-Rea(P)]*u(P)Im*E (u p ,p)*u(P), (22) 



where in the last line, we have used the effective Dirac equation, Eq. (21), satisfied by 
*u(P). uj p is the dispersion relation of the corresponding mode and is given, for example 
in Ref. [O]. Let us now define the spinor u(P) related to *u{P) by 



u{P) = [1 - Re a(P)] L/2 *u(P) , (23) 

thereby cancelling the gauge-parameter dependent factor in the normalization of *u(P). 
The spinor u(P) satisfies Eq. (j2"ID and is normalized according to Eq. (|TT]) . The explicit 
expression for the residues of the one particle contributions to the propagator is [|TT] 



Up-p 2 

Zp ~ 2m) ■ (24) 
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Eq. ( p2|) is now written in terms of u(P) as 
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— u(P) Im*Z (u p ,p)u{P), 



(25) 



which is manifestly gauge-parameter independent. Arguments to show that the expression 
for the damping rate in the Coulomb gauge coincides with that in a covariant gauge have 
been given elsewhere ||, therefore all of the above translates as well to the case of the 
Coulomb gauge. In the next section, we will write the explicit expression for the imaginary 
part of the gauge independent piece of the self-energy, *S . 

4 The reaction rate 

We now proceed to establish the relation between *£o and the total reaction rate T in 
the HTL approximation. In order to do this, we need an explicit expression for the 
imaginary part of *£o- The imaginary-time expression for the gauge-independent part of 
the self-energy is obtained from Eq. (|13|). To evaluate the sum, we use the identity (see 
for example Ref. [ |T2]j ) 



where the analytical continuation iuj — > po+ie has been performed. Here / and / represent 
the statistical distributions, whereas p and p are the spectral densities corresponding to the 
functions g and g, respectively. The spectral densities contain the discontinuities across 
the real energy axis. Their support depends on whether the momentum four-vector is 
inside or outside the light-cone. For time-like momenta, these spectral densities have 
support on the quasi-particle poles, whereas for space-like momenta they have support on 
the whole interval corresponding to the branch cut of the propagator. 

Adding up the two contributions in Eq. (|T3|) , Im*£ (u;p,p) can be expressed, by means 
of Eq. (P^|), as an integral involving products of spectral densities [|HJ, 



n 



x f^f ^f(k )f(p' )5(p -ko- p'oWM) , 

J-oo Z7T J-oo ZlT 



(26) 



lm*T, (u p ,p) 



x 




+ 7 m Pm^o, k)p(Pa,P- k)F?{P,K) 
+ Yp^(ko,k)a(p' ,p-k)F»(P,K) 
+ Fj t (P,K)p lw (k Q ,k)p(p' Q ,p- k)Y 
+ F»(P, K) Plu ,(ko, k)a(p' Q ,p- k)Y 
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+ F?(P, K) Pflv {k , k)p(p' , v- k)F?{P, K) 
+ F?(P, K) Plw (ko, k)a(p' ,p- k)F»(P, K) 
+ F$(P, tf)/V(*o, k)a(p' ,p- k)F»(P, K) 
+ if (P, K) PlMU (k , k)P(p' ,p- k)F»(P, K) 

+ G^{P,K)p» v {k Q ,k)} po=Wp , (27) 

where a g = g 2 /(4:7i) and the functions p^, p, a, $, F-f , F$ and G^ u are defined in the ap- 
pendix. Eq. ([27|) develops infrared divergences associated with the exchange of unscreened 
long wavelength magnetic gauge bosons. Nevertheless, it has been shown that the leading 
divergences can be resummed by a generalization of the Bloch-Nordsieck model at finite 
temperature |TJ]]. Eq. ( P7| ) is valid for any SU(N) gauge group, provided that the thermal 



masses include the corresponding group factors, and the generalization to the case with 
finite chemical potentials is also straightforward [fLlp . 

When computing the on-shell matrix element u(P)Im*T,o(ujp, p)u(P) using the medium 
spinors satisfying the normalization condition in Eq. (|IlD , the right-hand side of Eq. fl27f) 
can be identified as a sum of rates for processes involving thermalized fermions and gauge 
bosons. To this end, we note that the transverse and longitudinal projection operators 
correspond to the products of transverse (summed over components) and of longitudinal 
photon polarizations, respectively, and that 

(7° T 7 ■ P) = -ij: E uf(up,p)uf(^,p) , (28) 

U p Z p 3 

where the sum is over helicities of the (medium) fermion modes and the ± on the r.h.s. of 
Eq. (|28D refer to modes with either positive or negative helicity over chirality ratio. The 
functions F^, Frf and G^ v are effective vertices and the terms containing a or f3 represent 
products of amplitudes with an intermediate fermion propagator. Therefore, we write 

u{P)lm*i: {uj p ,p)u{P) = -UJ P T , (29) 

where T is the total reaction rate, which is the sum of decay and creation rates. Detailed 
balance ensures that the rate for the decay processes Vd, and the rate for the creation 
processes Tj, are related by T D = e^IY Eq. (jgD was first computed, for the case 
p = 0, in Refs. [12), Comparing Eqs. (|2~5|) and (^) we finally obtain the relation 



7 = T. (30) 

Three different kinds of products make up the total rate: the pole-pole, pole-cut and 
cut-cut contributions [|16j. The pole-pole terms represent processes involving two quasi- 
particles, a fermionic excitation with either positive or negative helicity over chirality ratio 
and a gauge boson excitation, either transverse or longitudinal, in addition to the original 
fermion mode. The pole-cut terms represent either scattering or annihilation processes 



9 



involving the exchange of a fermion or a gauge boson mode. Finally the cut-cut terms 
represent the same scattering or annihilation processes as the pole-cut terms together 
with the radiation of a soft gauge boson. The physical origin of the two latter kinds of 



processes is Landau damping |17]]. For specific values of the energy and momentum of 
the original fermion mode, some of these processes will be kinematically forbidden. It 
is easy to check that by restricting our attention to the pole-pole contributions in the 
first term of Eq. (|27D , one obtains the analogue to Eq. (35) in Ref. ||. We emphasize 
that, although our expression for T in Eq. fl2"9|) is formally the same as the one used in 
Ref. f2|, the original fermion mode is represented by the spinor u(P) corresponding to a 
quasiparticle propagating in the medium, and not by a free-particle spinor, from which it 
differs by the normalization. 



5 Hard axion production in a hot QED plasma 

As an application, consider the computation of the production rate of hard axions in a 
hot QED plasma through the Compton-like scattering diagram depicted in Fig. 2. This 
process is free from infrared divergences and thus it better illustrates the necessity of 
including the correct fermion wave function normalization in the calculation. 

The tree level calculation involves an integral over the three momentum transferred p' 
of the virtual electron. Since in the soft p' region HTL corrections to the electron propa- 
gator are not suppressed by powers of e, we need to use the effective electron propagator. 
Also, notice that since the axion momentum k is hard, then, for soft p', the momentum of 
the incoming electron p must also be hard. The contribution to the hard-axion production 
rate, T(k), from the soft p' region is most conveniently found from the imaginary part of 
the axion self energy depicted in Fig. 3. The heavy dot in the loop electron line represents 
the effective propagator. 




Figure 2: Axion photoproduction by Compton-like scattering. 
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In Euclidean space, the expression for the axion self-energy is written as 

n(*) = X 2 TJ2J 0-Tv[^ 5 S(P')S H (P' + K)^ 5 ] , (31) 

where A is the effective axion-electron coupling with units of inverse energy and Sh is the 
hard eletron propagator. Given that at least one of the momenta flowing into the vertices 
is hard, there is no need to consider HTL corrections to the vertices. T(k) is related to 
U(k) by 

1 Imll(fc) . . 

r « = MT^T (32 > 

where the factor (1 — e k °/ T )~ 1 selects the processes in which axions are produced. We now 
use the analog of Eq. fl26|) for the fermion-antifermion case to compute the imaginary part 
of the axion self-energy. The terms involved are 

lmTj2S(P')S H (P' + K) = ^(l-eWT) rM 

x f{p'o)f(-Po)S(ko + p'o - Po)p(p' )p(Po) , (33) 

where in the right-hand side we have performed the analytical continuation to Minkowsky 
space. Notice thatpl| since p, p are hard and pi — p 2 > 

P±(Po,P) = 2ixZ + (p)5(p t w+(p)) , (34) 

since for hard momentum, the modes with negative helicity over chirality ratio decouple. 
Keeping only the contribution from positive energy electrons, the hard axion production 
rate in the soft p' region is written as 

r -« - ^fwtflf /;^/«)«o +rt -.) 



x 27rZ + (p)5(p -u + (p))Tr[f(p(p' )fi(yo-j-p)}. (35) 




P 



Figure 3: Lowest order hard axion self-energy. The heavy dot in one of the internal 
electron lines represents the HTL effective electron propagator. 
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After integration over the direction of p 1 in Eq. ([35]), we obtain two kinematical restrictions 
for p' and p', namely p' 2 — p' 2 > and 2k > p' — p' Q} therefore, after some simplifications 

x {{3 + {p' ,p>)Tr[ft{ l0 - 7 • p')ft(jo -1-p)} 

+ P-(p' Q ,p')Tr[fi( l0 + j-p')fC( l0 -j-p)}} , (36) 



where we have introduced the intermediate scale q* with eT <^ q* r||18|. Carring out 
the calculation of the traces and using that for hard k 



/( Po k ) ~ e -(k +p' )/T + 1 ~ e - k /T + i ~ 1 ' ( 37 ) 



we get 



T soft(k) = dp' J_ pi dp>J(p' Q )Z + (p' + k)6(2k + p> -p')(p' 2 -p' 2 ) 

x {(p' -p' )(2k + p' +p')(3 + (p' ,p') + (p'+p' )(2k + p' -p')(3.(p' ,p')} (38) 

We notice that Eq. (|38"D involves a factor Z + (p' + k) which arises from cutting the electron 
line with hard momentum in the loop diagram of Fig. 3. This factor signals that the 
electron at hand corresponds to an excitation in the medium. 

We now ask whether the same expression for the hard axion production rate can be 
recovered from considering the self-energy of this hard electron in the medium. As we will 
show, the expression thus obtained will coincide with Eq. (|38|) provided that we include 
the proper normalization of the medium electron. 

Consider the diagram in Fig. 4 representing the self-energy of a hard electron in the 
medium. The dotted line represents the axion. As before, in the soft p' region, HTL 
corrections to the electron propagator are not suppressed by powers of e and have to be 




Figure 4: Lowest order hard electron self-energy. The heavy dot in the internal electron 
line represents the HTL effective electron propagator. 
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resummed into the effective electron propagator represented by the line with the heavy 
dot. In Eucledian space, the expression for the electron self-energy is 

S(P) = -A 2 T£ / ^(P+P>) l5 G(P + P')S(P')(P+P') l5 , (39) 

where G(P) is the bare axion propagator. Production of axions is equivalent to destruction 
of the original hard electrons. The rate of destruction of these electrons is given, according 
to Eqs. flTJ) and Q3IJ), by 

r(p) = -z +( p) ^ (1 - + ^™ , («) 

where the factor (1 + e~ PQ ' T )~ l selects the processes where the hard electron disappears 
from the medium and we have also summed over helicity states of the incoming electron. 
Notice that in Eq. ([TOD we have included explicitly the factor Z + (p), as corresponds to a 
rate that involves medium electrons. Using Eq. we have 

lmTj2G(P + P')S(n = n(e^ T + l) r^f(k )f(-p' ) 

n Joe Z7T Joo Z7T 

x 5( Po -k + p' )p(p' )27r5[k 2 - (p + f?) 2 } , (41) 

where we used that 21mG(P+P') = 2tt5[(P+P') 2 ] and have also performed the analytical 
continuation to Minkowsky space. From Eqs. (f|0j) and ([41]) we get after carrying out the 
angular integration 



(P) = Tfi - ^ [ 9 dp'p' f P dp' o e(2p + p' Q ~p')—^-e Po/T f(po+p' )f(~p' ) 

l07T JO J-p' p 

x Tr[( 7o - 7 -p)(P+nP{Po)(P +/")], (42) 



where the integration over p' is up to the intermediate scale q*. Notice that since po is 
hard and p' Q is soft, 

P Po/T po/T 

^ f ^ +p , ) ^ __ K _ = 1 (43) 

Thus, after computing the traces involved, we get 

* f q * iuu f p ' i„iom„ , ^ir ,,Z + {p)(p' 2 -p> 2 ) 



r«,/t(p) = 7^-4/ rfpV / dp e(2 P +p -p)f(-p y 

l07T 4 JO J-n' 



16vr 4 Jo ^ J-p> " u v * " u ™ p 2 

x {/5 + (p o ,p0(p / -Po)(2p-p , +Po)+/54Po ) P0(p' + Po)(2p + p' + Po)} (44) 
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In order to transform the rate of electron loss into the rate of axion production, we use 
that p = k — p' , as required by kinematics, and that k is hard. Therefore, 

1 1 1 



p 2 (k — p' ) 2 k 2 

Z+(p) = Z+(k-p' ) 
(2p + p' Tp') = (2k-p' Tp'). (45) 

We now introduce unity in the form 

1 = J d 4 k5{K -P-P'), (46) 

to write 

r*,/*(P) = / d " k5 ( K ~ p - P')^soft(k) (47) 
where T so f t (k) is given by 

A 2 ri* rp' i 2 

t( k ) = ^p^i J d P ]_ , d Pof(-Po) z +( k -Po) d ( 2k -Po ~P)(P ~Po 



soft\ 



(//;„./ ( -\>»)z- - \>^)v\Lh- - i> (} - j> it//" 

p 

x {( p '- p ' )(2k-p' -p')(3 + (p' ,p') + (p' + p ' )(2k-p' +p')(34p' ,p')} (.48) 



Eq. (|48|) coincides with Eq. d38|) as can be checked by means of the substitution p' Q — > — p' Q 
for which f3 + <-» j3_. 

It is now straightforward to compute the explicit expression for Y so j t {k). We use that 
in the soft p', p' Q region f(p' ) ~ 1/2 and that k ^> p , p' , therefore 

x {(p'-Po)P + (Po,P') + (p' + Po)P-(Po,P')} ■ (49) 



The integral is computed by using the sum rules satisfied by f3± [[□]]. The result is 

\ 2 Z+(k) \ q* 2 r°° 
F S oft(k) = —-j 1 — m) < — 5- - / dx \(x 3 - y+(x)x 2 - yl(x)x + y%(x)) Z+(x) 
2kn J [mj Jo LV ' 

+ (x 3 + y^(x)x 2 -y 2 _(x)x -y 3 _(x)^j Z_(x) } , (50) 

where we have defined the dimensionless function y±(x) = u±{x)/w,f and have extended 
the limit of integration in the second term to oo. The remaining integral in Eq. ([50]) is also 
a dimensionless quantity that can be computed numerically and its value is approximately 
1.3. 

We now proceed to compute the contribution to the hard axion production rate from 
the hard momentum transfer region. For this case, the tree level electron propagator 



14 



must be used. Notice that the tree level propagator can be recovered from the effective 
propagator by simply neglecting mj in the denominators of Eqs. (16T[) . In this manner, 
the approximate spectral densities become 

-SIS- <5i) 

Thus, the contribution to the axion production rate from the hard p' region is 

IW*) = ^^™ 2 f £ d P' C dp'J{p',) {P ' 2 ~f ) e{2k +p' - p') . (52) 

It is convenient to further decompose the hard region (p' > q*) into a (/) low (p' Q < q*) 
and a (II) high (p' > q*) frequency regions. In the low frequency region, we can still use 
the approximation f(p' ) ~ 1/2. The remaining integrals are readily performed and the 
result is 



^L(k) = ^^rn) (-g* 2 ) . (53) 



In the high frequency region, / provides the cutoff for the integrals. Recalling that q* < T 
we get 

r£&(*) = ^t^™? + 2kTln(l + e-^)} • (54) 

Adding Eqs. fl50D , ( [53"D and (|54l) , the dependence on the intermediate scale q* cancels, as 
it should, and the final expression for the hard axion production rate is written as 

= ^|±W m 4 (4 + ^ ln(l + e- k ' T ) - 0.65 1 , (55) 
kn J [mj mj J 

which is valid for k ~ T '. Inclusion of a possible electron chemical potential is straight- 
forward. Equation (|55|) represents the leading contribution to the photoproduction rate 
of axions in relativistic plasmas [|l^] at large k. Though this process is known to be sub- 



dominant as compared to nucleon bremsstrahlung of axions in this kind of plasmas [20 



a complete leading order treatment such as that of Eq. ( p5|) could have significant conse- 



quences for astrophysical processes in strongly magnetized relativistic plasmas f2T |. 

We should enphasize that in order for the rate obtained from the axion self-energy to 
formally coincide with the rate obtained from the electron self-energy, the latter has to 
be weighed with the factor Z + as corresponds to a medium electron mode. 

Finally, we note that for hard fcfTT| 
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thus, the error involved when ignoring the factor Z + is small. This is in contrast to the 
situation when k is soft for which 



Z ± (A0«i±-^-. (57) 
2 6m f 



6 Discussion 



In this paper, we have argued, in the realm of hot gauge theories, the equality of the 
damping rate of a fermion mode propagating in a medium and its total reaction rate, 
to leading order, provided the latter is computed using the wave function of the mode 
in the medium. This result is based on the explicit gauge independence of the on-shell 
matrix element of the absorptive part of the fermion self-energy, according to the effective 
perturbative expansion of Braaten and Pisarski. It also unifies two apparently distinct 
concepts, both related to Im*£ . In addition, Eq. (|30|) also sheds light on the distinct 
interpretation of a quantity like T in the work of Wei don [fj] and in ours. In effect, 
according to Eq. (]TT|), the spinors u(P) are normalized to Z p , related to the probability 



of finding a fermion mode with energy u p and momentum p pertaining to the medium. 
On the other hand, if the spinor u(P) was taken as free, it would represent a mode in 
vacuum. Thus, if we take spinors from the medium, T can be interpreted as the rate to 
approach equilibrium for an excitation in the medium — slightly out of equilibrium - 
which, according to Eq. fl30|), is equivalent to the rate at which the excitation is damped. 
On the other hand, if the spinors are from vacuum, the corresponding reaction rate can be 
interpreted as the rate to approach equilibrium for a test particle, taken from the vacuum 
and dropped into the medium. 

We have used the equivalence between damping and reaction rates of fermions in a 
medium to compute the hard axion production rate in a hot QED plasma by Compton- 
like scattering processes. We have shown that the rate obtained from the axion self- 
energy coincides with the rate obtained from the electron self-energy only when the proper 
medium wave function is used in this latter computation. 

Another interesting point, that we do not touch upon here, has to do with whether 
Eq. ([ft]) is a general result, independent of perturbation theory, which is related to the 
way the thermalization of gauge degrees of freedom is treated P2"| . The fact that at least 
to lowest, but complete order, the expression for the fermion damping rate in hot gauge 
theories is identical in the Coulomb or in covariant gauges, hints towards the validity of 
the result to all orders. 
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Appendix 

Here, we list all of the functions introduced in section IV. Capital letters are used through- 
out to denote momentum four- vectors, all of them in Minkowski space and = (— 1, q). 
We start with the functions Fi, F 2 and G^, related to the explicit expressions for the 
HTL effective vertices 17, and 17 
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F?(P,K) 
F£{P,K) 
G^ V (P,K) 



—m. 



—m. 



-m. 



cm Q^Q 



V 



1 



4vr[p.Q] \[(P-K)-Q], 



4vr [P ■ Q] 

[P ■ Q] 2 



dtt 



5{{P-K)-Q) 
5{{P-K)-Q). 



(5f 



V denotes the Cauchy principal value These functions depend on the quantity 

A = P 2 (P - Kf - (P ■ (P - K)f and vanish for A > 0. 
Next, let us define the spectral densities, p^ v , p, a and /3, 



P~(Po,P) 
f3(Po,P) 



k' 2 



k 2 



;p L {k Q ,k)Ph, + pr{ko,k)Pl 



P+(Po,p)(lo ~ 1 ■ P) + P-(po,p)(lo + 1 ■ P) 
«+ (Po, P) (7o - 7 " P) + «- (Po, P) (7o + 7 ■ P) 
(3+ (po, P) (7o - 7 • P) + P- (Po, P) (lo + l-p): 



(59) 



where P L and P T are the longitudinal and transverse projection operators in four dimen- 
sions. The quantities Pl,t, P± are the standard HTL spectral densities corresponding to 
the gauge boson and fermion propagator, respectively |TT|| . We have also defined 



a±(Po,p) 



(27T) 



PoTp-^(^ln 



PO+P 




po-p 





PO+P 

Po-P 



± 1 



+ 



™7 A T PO 

2p V P 



n 2 ' 



(60) 



P±(,Po,p) 



^ 3 ?(™)% 2 -Po) 



PQ+P 
po-p 



± 1 



+ 



2p V + P 



(61) 
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